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Abstract 

In this paper, we present a coherence protection method based upon 
a multidimensional generalization of the Quantum Zeno Effect, as well 
as ideas from the coding theory. The non-holonomic control technique is 
employed as a physical tool which allows its effective implementation. The 
two limiting cases of small and large quantum systems are considered. 

1 Introduction 

The uncontrollable interaction of a quantum system with its environment is 
responsible for "quantum errors" which lead to a partial or complete loss of 
the information initially stored in its quantum state. After Shor's demonstra- 
tion ^ that error-correcting schemes exist in quantum computation, a general 
framework of error-correction has been built upon the formalism of quantum op- 
erations. The main contributions concern quantum codes, and particularly the 
class of stabilizer codes ; other strategies developed suggest the use of " noiseless 
quantum codes" or " decoherence-free subspaces". All these methods usually 
demand that errors act independently on different qubits (the independent er- 
ror model), and make use of the symmetry properties associated with these 
requirements. This implies that the set of errors to be corrected is restricted 
to a special subgroup, called the Clifford group. In this paper, we present a 
protection method which does not make so drastic assumptions. 

For low dimensional systems, one can take advantage of the Quantum Zeno 
Effect allied with basic ideas of coding theory in order to protect the infor- 
mation encoded on a subspace of the total state space. To be more explicit, 
one frequently repeats a three step sequence comprising coding, decoding and 
projection, which prevents errors from developing in the system : coding and 
decoding consist in the application of a unitary matrix and its inverse, which 
can be achieved through non-holonomic control, and act in such a way that 
erroneous infinitesimal components are orthogonal to the information subspace 
; projection is performed through an irreversible process such as spontaneous 
emission which clears the unwanted orthogonal increments out. 

Promising for small quantum systems, this method becomes however expo- 
nentially complex when the number of the qubits involved increases. For large 
quantum systems, we suggest to employ random coding to reduce the influence 
of errors consisting in binary interactions. In this context, high dimensionality 
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does not appear as an impediment, but rather as an advantage, since it 'dilutes' 
the influence of the errors. 

This paper is organized as follows. In the second section, we show that the 
Quantum Zeno Effect allied with basic ideas of coding theory allows one to pro- 
tect the quantum information contained in low-dimensional quantum systems. 
In the third section, we present the random coding technique which can protect 
the information stored in large systems against the errors resulting from binary 
interactions. 

2 Coherence Protection in low-dimensional sys- 
tems : Quantum Zeno Effect and Non-Holonomic 
Control 

The Quantum Zeno Effect (QZE) [3131 appears in a system which is frequently 
measured in its (necessarily known) initial state: if the time interval between 
two projective measurements is small enough, the evolution of the system is 
nearly "frozen" ; in other words the one-dimensional subspace spanned by the 
(necessarily known) initial state is protected against the influence of the nat- 
ural Hamiltonian of the system. We suggest to generalize this effect in order 
to protect any (unknown) vector of a given multidimensional subspace of the 
whole Hilbert space. To this end, we propose an information protection scheme, 
described in the first paragraph of this section as well as the algorithmic 
tools which allow its implementation, and which are presented in the second 
paragraph. 

2.1 Multidimensional Zeno Effect and Coherence Protec- 



In this paragraph, we shall first present the multidimensional QZE which allows 
us to protect an arbitrary subspace of the Hilbert space against the action 
of a set of given interaction Hamiltonians. Then, we shall take advantage of 
this phenomenon to protect an information-carrying subsystem of a compound 
quantum system from the influence of some uncontrolled error-inducing external 



Consider a quantum system S, whose A'^-dimensional Hilbert space is de- 
noted by Ti. and whose time-dependent Hamiltonian has the form 



tion 



fields. 



M 




(1) 



m — 1 




are M given independent Hermitian matrices on TC and 



{fmir}} 



m=l,...,M 



are M unknown functions of time. The Hamiltonian H{t) 
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accounts for the errors we want to get rid of. Note that the unperturbed part 
of the Hamiltonian is assumed to be zero (or proportional to the identity so 
that one can set it to zero). The standard QZE allows us to nearly "freeze" the 
evolution of the system by measuring it frequently enough in its (known) initial 
state ; in other words, through this effect we can protect the one-dimensional 
subspace spanned by the initial state of the system from the influence of the 
error-inducing Hamiltonian In what follows, we generalize this effect so as 
to protect an arbitrary multidimensional subspace C from H{t). 

Any vector of C evolves according to U{t, to) — T |exp —i J^^ H{T)dT | , 

where T denotes time-ordering, and where we set h—\. For the QZE to hold, 
we shall only consider evolution in short time periods, the duration tz of which is 
so short that the corresponding action of the M components of the Hamiltonian 



CP) is small, i.e. Em //^"^ fm{T)dT 



<C 1. We can thus expand 



= 1 



Uit + Tz,t)^ U,nf -I-^y^[ fm{T)dT E^- (2) 



After a Zeno interval r^, the initial state \ip) is thus transformed into lip^) = 

\tp) + \Stpe) where IStpe) ^ -«Z]m=l ^mEm \lp) with Ejn = (/ fm{T)dT) . 

Let us assume that we are physically able to perform the measurement- 
induced projection onto C in the system S (see below the discussion of such 
projections for compound systems comprising an information subsystem and an 
ancilla). If we straightforwardly apply the standard QZE procedure by merely 
projecting the state vector \ipf,) , resulting from the infinitesimal evolution of the 
initial state lip) , onto C, we get the vector \tpp), which, a priori, differs from 
(see FigOJi), since, usually, the vectors Em IV") and thus the increment vector 
\Sipe) itself are not orthogonal to C. It is thus clear that we have to adapt the 
standard Zeno strategy. ^ 

To this end, we assume a unitary matrix C acting on 7i, which we call 
the coding matrix, such that the Hermitian operators Em > act or- 

L J m=l,...,M 

thogonally on the subspace C = CC, called the code space. Let us denote by 
J > 1 the dimension of C and by {|7i)}j^]^ / one of its orthonormal bases ; 

< |7i) = C \ > will denote one of the orthonormal bases of C, the state 

vectors being called the codewords. For any pair (1%) , \jt)) of codewords 

and any operator Em € \ Em r we have, by the definitions of C and C 

I J m=l,...,M 

(7t|7s) = Sst (orthonormality conditions) (3) 
( 7t Em 7s ) = (orthogonality of the errors) (4) 



Equivalently, for any pair (jV') , |x)) of vectors of C and for any operator Em € 
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L J m=l, 



X 



C^E^C 



0. 



(5) 



In particular, for any pair {\js) , \lt)) of basis vectors of C and for any operator 

Em G \ Em )■ 

I J m=l. 



It 



EmC 



Is 



0. 



(6) 



If we apply the coding matrix C to the initial state vector , before exposing it 
to the action of the Haniiltonian J^l, we obtain the new vector V'") = C IV") € C 

(Fignj3l,2) which is transformed after a Zeno interval tz into 



, where 



Decoding 



-06/ yields the vector \ Ae) ~ C 



^ ~» E™=1 £mEmC IV) (FigHtS) 

= \A + \SAe) where \6Ae) - 



i J2m=i £mC^E.mC {ip) ■ From Eq.|(SJ) it can be seen that for any vector \x) G C, 



EZ=i (xl C^EmC IV) = which means that ISip'e) is orthogo- 



nal to C (Fig^4). A measurement-induced projection onto C finally recovers 
the initial vector \ip) with a probability very close to 1 (the error probability 
is proportional to rf ). If the coding-decoding-projection sequence is frequently 
repeated, any vector {ip) of the subspace C can thus be protected from the 
Hamiltonian for as long as needed. 

The multidimensional generalization of the QZE we have just presented al- 
lows one to protect any subspace C of a Hilbert space Ti. against Hamiltonians 
of the form and is thus very useful in the context of information protection 
as we shall see in the following. 



Indeed, let us consider an information system X of Hilbert space Tij and 
dimensionality /. This system is subject to a set of M error-inducing Hamil- 
tonians \ Em \ which, for instance, represent the interactions of the 

L J m=l,...,M 

system with M uncontrolled external classical fields fm{t)- we want to get rid 
of this external influence which is likely to result in the loss of the information 
stored in the initial state vector |V/) = X]i=i Wi), where {|'^i)}i=i / denotes 
an orthonormal basis of Hj. To this end, we shall use the multidimensional 
Zeno Effect. As the multidimensional QZE can only protect a subspace of the 
whole Hilbert space, we first have to add an A-dimensional auxiliary system A 
(called ancilla) to our system T, so that the information is transferred from Hi 
into an /-dimensional subspace C of the (N = I x A)-dimensional Hilbert space 
H = 1-Li®1-La of the compound system S = X ® A. Furthermore, we shall sup- 
pose that all the state vectors of the different Hilbert spaces H/, H-a and hence 
H. are degenerate in energy so that the unperturbed part Hq of the Hamiltonian 
can be set to zero as in the first part of this section: the subspace C and the 
information it carries can thus be protected through the multidimensional QZE. 
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a) Without coding-decoding 




b ) With coding-decoding 




Figure 1: Multidimensional QZE: a) a simple projection fails to recover the 
initial vector, b) the sequence coding-decoding-projection protects the initial 
vector. 
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Note that A and X need not be "physically separate" systems, but only have to 
possess independent Hilbert spaces Ha and Hj. 

Let us now return to our problem and first consider the simple case when the 
ancilla is initially in the pure state The information initially stored by \ipi) G 
Hi is transferred into the factorized state ji/') = ® \a) — = 



Equiva- 



ELi c^ hi) of C = Hi®Span [\a)] = Span \{h,) = (g) |a)};^i 
lently, the initial density matrix of the compound system S is p = {lipi) <E) 
(la) which is transformed after the coding step into p = C^pC ; at the end 
of the action of the errors it is transformed into p^ ~ Ul^^C^pCUmf ; finally, 

after decoding, it takes the form pe — CUj^fC^pCUinfC^- In this setting, the 
projection onto C can be simply achieved by measuring the ancilla in its initial 
state \a). As tz is very short, the state of the ancilla evolves just a little within 
a Zeno interval : the probability of detecting it in its initial state \a), and thus 
of projecting the state of the compound system onto C is thus very close to 1. 
After projection, we trace out the ancilla to obtain the final reduced density 
matrix p'^ — {a\ CUj^^C^pCUinfC^ \a) for the information system X; in the 
same way, one can calculate the initial reduced density matrix is pi — \ipi) {ipi \ ■ 
The variation Spi = p'j — pi of the information-space density matrix during the 
whole process can then be expressed as the commutator 



Spi = -i 



M . 

\a) ,pi 



from which we infer that pj satisfies the equation — 



i,Pi 



where 



dt 

he — X)m=i frn E^C \a) IS an effective Hamiltonian which is determined 

by the error-inducing Hamiltonians transformed by the coding and decoding 
and projected onto the initial state of the ancilla. From Eq.© one can see 
that he = Q and hence pj remains constant in time: as long as we repeat the 
coding-decoding-ancilla resetting sequence, the information initially stored in I 
is protected. 

It is not always feasible to directly measure the ancilla independently from 
the information system ; in other words, it is sometimes impossible to perform 
a projection onto disentangled subspaces of H of the form Hi®Span [|a)] : in 
some cases, one can only project onto entangled subspaces of the total Hilbert 
space H. In such a case the information initially stored in the vector = 
Si=i hi) ^ transferred into an entangled state of I and A of the 

form h) = where the / vectors |7i) (z = 1,...,/) which form an 

orthonormal basis of the information-carrying subspace C, are not factorized 
as earlier but are in general entangled states. Nevertheless the same method 
as before can be used in that case to protect information, albeit in a different 
subspace C. 

To conclude this description of our method, let us now return to conditions 
^ and ^ imposed on the codewords {[7^) , z = 1, /} and make two remarks 
about them: 
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A. We can establish a useful relation between the dimension A of the ancilla 
and the number M of correctable error Hamiltonians. The set of the / codewords 
can indeed be seen as a collection of 2/ x iV = 21^ A real numbers on which 
2/2 + 2MP = 2/2(1 + M) constraints, directly derived from Eas. pH) . are 
imposed. As the number of free parameters must be larger than the number of 
constraints, we necessarily have 2PA > 2/^(1 + M), or equivalently 

A~1>M. (7) 

This condition, called the " Hamming bound" , gives an upper-bound on the 
number of independent error-inducing Hamiltonians that our method can cor- 
rect simultaneously. 

B. We may compare our correctability conditions Q with the more general 
conditions of standard quantum error-correction jJl |H] 

V(|7s),|7t» GC2, v(Efe,E,) e{E, ({^„})}, 

7s) = "fe; (7t|7s) (8) 



It 



which ensure the existence of a code space that is completely protected against 
the error-inducing Hamiltonians Em- Here an are complex numbers, and the set 
|i?m| of Hcrmitian operators Em generates a group G ({-^™}) '^^ possible 
error- induced evolutions By |Ej (|^'"}) } denote a complete basis set 

of operators which spans the space of evolution operators U and allows one to 
represent any U as a linear combination of the basis operators Ej. In addition 
to all the Em, the variety of all linear combinations of Ej includes also many 
other operators given by commutators of all orders in Em entering the expansion 
of U for long times. The condition ||SJ) is therefore much more restrictive than 

Eq.|01). Moreover, even for two generic matrices Em, the basis spans the 

entire Hilbert space Ti., yielding C — Only if the set |i?,n| belongs to an 
extraspecial algebra restricting the error evolution operators C/ to a subgroup 
Q ({-^m}) i'H) of the full unitary group in H, a non-trivial code space C 

may exist. The Zeno effect is the only way to suppress loss of coherence if it is 
not the case. 



2.2 The code space and the coding matrix 

It is sometimes possible to build the code space C explicitly from physical con- 
siderations. However, in general, we need an algorithm to calculate the code 
basis {|7i)}j=i / or, equivalently, the coding matrix C. In this paragraph, we 
shall first describe this algorithm, then, we shall show that the non-holonomic 
control technique can be employed to implement the coding matrix physically. 
We will also provide an algorithm which achieves the appropriate control. 
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Let us first make a remark which will be useful. Consider a vector |C) of 
some Hilbert space and a matrix E on this space. From the vector |C) we want to 



calculate a vector 

and the function /^(A) = 
minimal for A = 0: indeed 



C ) such that ( C 



AE 



C) = 0. If (C 



C) = 0, then |C) 



depending on the c- number A, is 



AE 



CIC) + A(C 



= 1 + |A|'(C 



EtE 



A* (C 



Et 



C) + |A|VC 



EtE 



and as (C 



C 



EtE 



C 



> 0, /f;(A) is minimal for |A| = 0, that is A = 0. But, if 
C) 7^ 0, we can apply the following iterative method: we minimize /c(A) 



with respect to A, then we set |C') = |C) + |C) and take 



as our new 



C) = 0. 



|C) ; repeating this sequence finally leads Cy, such that i^C 

Let us now return to our problem and show how to use the previous re- 
mark. We want to find / vectors which meet the conditions Q and Q ; 
equivalently, we look for an orthonormal basis in which all the matrices Efc have 
their I x I upper left blocks equal to zero. To solve this problem, we propose 
to transform our initial problem in such a way that it can be dealt with by the 
iterative algorithm presented in the previous paragraph. Let us combine the / 
vectors |7i) into a {N x I) " supervector" 




Then let us build E = (-^^^ + M^^-^^ different (iV x /) x (A^ x /)-dimensional 

super-matrices E^ in the following way: we consider them as made of blocks 
of dimension N x N and we successively fill each of these blocks with the dif- 
ferent Hamiltonians Em or the identity matrix / or 0. To be more explicit, 
the first ^'"''2^'' matrices are built by simply placing the N x N identity matrix 
in each of the ^^'^2^'^ blocks situated above the diagonal. In the last ^^^^'^^') 
ones, the M operators £',„ are successively placed in each of the ^^^^^^ blocks 
on and above the diagonal. One can thus reformulate the conditions (j2Jl as 



follows: for 1 < fc < 



C 



C ) = 0. This form docs not take the 



normalization condition into account, which will be imposed in a different man- 
ner. Similarly, the conditions Q are translated into the following form: for 



1 < fc < 



A//(/+l) 
2 ' 



0. This new problem can be 



2 I J- ^ ^ 2 ' 
handled by the same kind of iterative algorithm as in our preliminary remark. 
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First, we randomly pick a supervector |Co) which will be the starting point of 
the first step: we normalize this vector by imposing to each of its / components 
to have norm = j. If one of the components of |Co) is non normalizable, that 
is equals zero, we pick up a new random supervector |Co) as a starting point. 

Then, we minimize Fq, (^xf\ \^°\ \f) - Ef=i || |Co) + A^"^ |Co 

with respect to the E c-numbers xf \ and we calculate |ACo) = J2k ^'f'^k |Co) 
and |Co) = |Co) + ^ |ACo). We normalize jCg) by requiring each of its / com- 
ponents to have the norm = j , and take the result of this operation as our new 
starting point |Ci). If one of the components of \Cq) is non normalizable, that 
is equals zero, we pick up a new random supervector |Co) as a starting point. 

We repeat this sequence of operations as long as needed and obtain the de- 
sired vector C ) asymptotically. Practically, as our algorithm converges quickly, 



the number of iterations needed is small. 



The coding matrix C is a complex unitary operator on the Hilbert space 
of the compound system S = T ® A. We have just shown how to calculate the 
codewords, which actually form the first / columns of C, but one can wonder how 
to implement it physically. This question can be solved by the non-holonomic 
control technique. 

Indeed, we can directly apply the results of the first of our articles to our cod- 
ing problem in the following way: first, we find the codewords , « = 1, /} 
by the iterative algorithm we have previously presented, then we complete the 
set of / vectors {|7j) ,i — 1, 1} with {N — I) vectors {|7j) , j = I + 1, N} 
to form an orthonormal basis of 7i, we build the coding matrix by taking the 
vectors {|7i) ,i = 1, ■■■,N} as columns of C, and finally we calculate the N'^ 
appropriate timings {r^} such that 

U (ti, ■■■,tn2) = exp {-iHaTm^ ■ ■ ■ exp (^-iHhT^ = C 

through the complete control algorithm we have previously presented (we sup- 
pose we have two distinct perturbations Pa and Ph such that the system is 
completely controllable). Note that we assume Hq = 0, hence Ha — Pa and 

Hb = Pb- 

Actually, this straightforward procedure provides a lot of useless work. In- 
deed, most of the information contained in the coding matrix is irrelevant and 
the A'^^ real parameters of C do not all have to be controlled exactly: the num- 
ber nc of necessary control parameters {r^} is actually much less than A^^, as 
we shall see now. 

The coding matrix is characterized by the relations The problem of 
control thus reduces to finding nc timings r^, forming the time- vector r = 
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such that the non-holonomic evolution matrix 



U ( T*) — exp I 



{^-iHaTnc^ . . .exp 



checks (jnj- The number nc of control parameters must exceed the number of 
independent constraints which is clearly ~ MP , that is nc > MP . Thus the 
number of necessary control parameters appears to be much smaller than iV^. 
So we need a new algorithm which achieves a partial and less expensive control 
of the evolution operator of the system. 

The algorithm we shall use to calculate the appropriate control timings Tj 
mixes the iterative algorithm presented at the beginning of this paragraph and 
the non-holonomic control technique. If we introduce the (TV x /) x [N x /)- 
dimensional block-diagonal matrix 





V 








and the (TV x /)-dimensional supervector 



|C) = 



/ l7i) 



V l7/> 



composed of the coordinates of the / basis vectors of C, we can set our problem 
of control into the following equivalent form: we look for a time- vector ~t such 
that 

Vfc, (c|ut(y)EfeU(y)|c) =0 (9) 
where the matrices \^k\ denote E different matrices of dimension [N x /) x 

L J k=l,...,E 

(TV x I) which have been introduced in the beginning of this section. The idea 
of our algorithm is to take the super vector |Co) = U (Vq) |C), where is a 
random time- vector, as the starting point for an elementary step of the iterative 

algorithm and look for the small time increment cLtq such that U a + ^'''"0^ |C) 
follows the direction provided by the result |Co) + I^Cq) of the iterative algo- 
rithm. The repetition of this sequence finally yields ~t — ~t q -\- cItq + dri + ... 
which meets Eq.(|5J). 

Let us now describe the algorithm in more detail. First, we randomly pick 
a set of timings ro,i in a " realistic range" , dictated by the system under consid- 
eration: in particular, control-pulse timings have to be much shorter than the 
typical lifetime of the system but much longer than the typical response delay 



11 



required by the experiment. Then we calculate |ACo) = J2k ^k^k |Co) by mini- 

, A^'' ) as in the algorithm presented 



mizing the same function Fqo ^A^' 



(0) ;^(0) 



at the beginning of this section. At that point, we look for the small increment 
dTQ of the time-vector ~tq such that 



Vfc, ( C 



1^ (^o) -d^o 



EfeU(yo) + ut(yo)Efc 



Co + iACo Efc Co + iACo)- (Co 




(Co + iACo|Co + iACo) 



(10) 



It should be noticed that we do not consider the error super-matrices E^, cor- 
responding to orthonormality conditions: in other words, we just take matrices 
|Efe| ,,,,, ,,,, into account. Thus we deal with ^^^(^+^) 



com- 



plex equations. This set of equations can be reduced to the real linear system 



Sir„)-dTo = Wi\ACo)) 



(11) 



where S'(Vo) and W(|ACo)) are respectively an MP x nc real matrix and 
a M/^-dimensional real vector. We obtained Ea. HllI) by splitting the set of 
^^^^"""^^ complex equations l(TU|l into two sets of ^'^^(_^+^) real equations, and 
rejecting those which are trivial (0 = 0) or redundant. Though straightforward, 
the explicit expressions of the different elements of S and W involve many 
indices and are so unpleasant that we prefer not to reproduce them here. 

The linear system we have just found is, a priori, rectangular {MP x nc), 
but actually the number nc has not been fixed yet. Previously, we stated that 
no > MP: we could be tempted to set nc — MP so as to obtain a square 
system, easily solvable by standard techniques of linear algebra. Yet we will 
proceed in a slightly different way. We set nc = MP + 5n > MP, where 
Sn is an integer of order 1, then we randomly pick MP timings ti among the 
nc which will be considered as free parameters, whereas the other Sn ones will 
be regarded as frozen. The new version of Eas ljllll is now clearly a square 

system, which yields the M/^-dimensional increment drg, corresponding to the 
MP free varying timings, which we complete with 5n zeros, corresponding 
to the frozen timings, into a np-dimensional vector drg. Then we set V" 
= ~tq+ a dTQ where a is a convergence coefficient and calculate the test function 

gU) - Efc (c ut(y)EfeU(y)' 



in r = T I for different values of 

a £ [0, 1]. If we find an ai such that G (V") < G (Vq), we take T*i = ~t'"^ 
as our new time-vector, and keep the same free-varying timings. If wc cannot 
find such an ai, this means we are situated in a local minimum of G ; then 
we set ~Ti = T*o and pick a new set of free varying parameters. This rotation 
procedure among control parameters allows us to avoid possible local minima 
of the test function G wc want to cancel. 
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We repeat this sequence of operations as long as needed and obtain the 
desired vector 't asymptotically. Practically, as our algorithm converges quickly, 
the number of iterations needed is small. 

We have not said anything about decoding so far. If the signs of the two 
Hamiltonians Ha = Pa and Hf, = Pj, can be reversed b^ altering the control field 
parameters, decoding amounts to reversing Pa and Pb and applying the same 
control timing sequence backwards. Otherwise, one must use the general non- 
holonomic control technique, involving iV^ control parameters, to find timings 
which realize C^^. 

3 Zeno Coherence Protection by Random Cod- 
ing 

The protection method we presented in the previous section seems promising 
for relatively low-dimensional systems. However, for large systems, it is likely 
to lead to very heavy computations and long control sequences. To deal with 
such systems, we therefore propose to employ an approach inspired by classi- 
cal random coding jTj : in this method, linear codes [n, k, d] are produced, in 
which fc-bit words are encoded as randomly chosen n-bit sequences. The mini- 
mal Hamming distance d between any two codewords approaches the Hamming 
bound d<n — kasn^oo. In this section, we show how to extend the idea of 
random coding to the quantum case. 

Strong mixing or entanglement occur in the phase or Hilbert spaces of com- 
plex classical or quantum systems, respectively, and can, in principle, be used 
for random coding. However, in practice, in the classical case, the dynamics of 
such systems is not reversible, which makes subsequent decoding hardly pos- 
sible. By contrast, the dynamics of multi-dimensional quantum systems can 
be reversed, when the underlying physical mechanism is simple enough : the 
spin-echo phenomenon is a typical example of this topic. High dimensionality 
of simple quantum systems is responsible for the massive parallel computing 
capacity of quantum computers. Therefore, we have to find an operation which 
produces strong mixing in the multidimensional Hilbert space, and which can 
be inverted in a simple way : the non-holonomic control suits perfectly this pur- 
pose. The essential requirement for the protection scheme we propose to apply 
is that error-inducing interactions are relatively simple, resulting, for instance, 
either from a binary qubit interaction or, generally speaking, from a few-particle 
coupling. 

To combine strong mixing with irreversibility, we assume that we have a 
quantum system with a large number of separate energy levels and with two sim- 
ple interactions which satisfy the bracket generation condition and can therefore 
be employed for the non-holonomic control. In such a system, one can encode 
quantum data into strongly mixed states by straightforwardly applying a uni- 
tary transformation C, the decoding procedure being achieved by the inverse 
transformation C^^ . Encoding the data into many levels allows us to strongly 
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reduce the error-rate, as will be shown. In turn, by applying the Zeno effect 
12 |3j we can restore the slightly corrupted data back to its original value with 
high probability. 

To encode the data in a high dimensional Hilbert space of n qubits, we intro- 
duce n — k ancillary qubits in addition to the required number k of data qubits. 
This results in an increase in the number of possible errors, which depends 
polynomially on n, but the error rate decreases at will because the infinitesi- 
mal errors are semi-orthogonal to the encoded data to a degree exponential in 
n — k. The degree of semi-orthogonality reflects the error correction efficiency 
of the coding. Efficiency requires a precise and careful choice of the code for 
coding in minimal dimensions, but in high dimensions it is naturally achieved 
by random coding. In mathematical terms, the method relies on the fact, that 
in a multidimensional space, a pair of randomly chosen vectors are almost or- 
thogonal with high probability. In physical terms, random coding is equivalent 
to strong mixing, or full population of all energy levels, which can be reached 
by the non-holonomic control with a number of interaction switchings depend- 
ing only polynomially on n. Thus the random coding approach of the present 
section complements our earlier non-holonomic Zeno coherence loss suppression 
scheme, which requires exponential effort to find and achieve an exact code, and 
is thus efficient only for low-dimensional systems. 

The essence of random coding can be elucidated as follows. Consider an 
n-qubit system, comprising a /c-qubit information carrying subsystem and an 
(n — fc)-qubit ancilla. In the 2"-dimensional Hilbert space of the system, the 
error-inducing Hamiltonians Em corresponding to a few-particle interaction are 
represented by sparse matrices in the computational basis, which is composed 
of all the possible tensor products of individual qubits eigenstates. In this 
basis, the number of non-zero matrix elements is indeed polynomial in n : for 
instance, for binary interactions, this number scales as . The coding-decoding 
transformation EmC, where C stands for a generic unitary matrix, 'smoothes' 
all the matrix elements by mixing them : finally, all these elements are of 
the same order of magnitude, which is, up to a polynomial factor, 2"-times 
smaller than the typical value of non-zero matrix elements in the computational 
basis before the coding-decoding sequence. Error matrices elements are thus 
exponentially reduced : the error suppression condition Eq. @ for the projection 
of these matrices onto the subspace of the initial state of the ancilla is not 
fulfilled exactly any longer : the projection differs from zero, but its norm 

^7f ErnC j-emaius small, and decreases exponentially with the 

size of the ancilla. The error-accumulation rate is thus inhibited by a factor of 
the order of the ancilla Hilbert space dimension. 

Note, that this mechanism is efficient only when the generic coding matrix 
C can be achieved by a small number of switchings, such that the coding pro- 
cedure does not take exponentially long time. Fortunately, the coding matrix 
takes a generic form after a relatively small number of switchings, which scales 
linearly with the number of qubits n (see |H]). Moreover, if the signs of the 
two interactions employed for the non-holonomic control can be inverted, the 
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decoding operation can be performed at the same level of complexity as the 
coding procedure, by straightforwardly changing the signs of the interactions 
and inverting the timing sequence in which these interactions are applied. The 
main restriction to practical implementation of the random coding protection 
scheme arises from that one necessarily has to remain in the Zeno regime: the 
measurement time, which does not depend on the size of the system, has to 
be much shorter than the coherence loss timescale, which decreases, although 
polynomially, with the size of the system. 

4 Conclusion 

The non-holonomic control allied with the Quantum Zeno Effect can be em- 
ployed to overcome the influence of the environment on the quantum system 
considered. On the one hand, in the case of low-dimensional systems, we 
showed that quantum information can be protected by frequently repeating the 
cycle coding-infinitesimal errors-decoding-projection : coding and decoding cor- 
respond to a imitary transformation of the Hilbert space and its inverse, respec- 
tively, which are determined in such a way that the projection onto the initial 
information carrying subspace of the state resulting from coding-infinitesimal 
errors-decoding yields the initial state vector. All the needed algorithmic tools 
have been presented. On the other hand, for high-dimensional systems, one 
can adapt the classical idea of random coding to the quantum case : the basic 
principle is to use non-holonomic control to impose generic and easily reversible 
unitary evolutions to the system in order to encode/decode the information ; 
this procedure "dilutes" the influence of the errors in the large Hilbert space 
and then decreases their influence. 
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